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O ! 

. Assuming that the order parameter corresponds to an equal spin triplet pairing symmetry state, 

£\j ■ we calculate the angular dependence of the upper critical field in triplet quasi-one-dimensional 

j i ' superconductors at high magnetic fields applied along the yz plane. 

o : 

The upper critical field of quasi-one-dimensional superconductors for perfectly aligned magnetic fields along the y 
(&') axis has been calculated for both singlet and triplet states . More recently new experiments performed in these 
systems lead to the observation of unusual superconductivity at high magnetic fields ^,^] and created new excitement 
in the area of organic superconductivity In quasi-onc-dimensional superconductors (of the Bechgaard salts 

family with chemical formula (TMTSF^X, where X = C104,PF6, ...) the upper critical field exceeds substantially 
the Pauli paramagnetic limit This indicates that high magnetic field superconductivity in these systems is 
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most likely triplet. The exact symmetry of the triplet phase and the mechanism of superconductivity are not yet 
| , known, however an early suggestion of triplet superconductivity (at zero magnetic field) was given by Abrikosov |Tc| l 
Oh' in order to explain the suppression of superconductivity in (TMTSF) 2 X in the presence of disorder (non-magnetic 
2 ; impurities) 

Very recently, it has been argued by Lebed, Machida, and Ozaki (LMO) that the spin-orbit coupling must be 
strong in order to explain the observed experimental upper critical fields and the absence of the Knight shift for fields 
parallel to the y (b') axis, [[uj However, we find this possibility hard to be realized in the (TMTSF) 2 X family since any 
reasonable estimate of the atomic spin-orbit coupling leads to a very small value (in addition, the Bechgaard salts are 
very low disorder systems, with very long mean free paths). Estimates of the value of spin-orbit coupling are several 
orders of magnitude smaller (jl5| than the values required to fit the critical temperature of (TMTSF) 2 PF6 ]l3| . |l6[ , 

■ for instance, at low magnetic fields. That the atomic spin-orbit coupling should be small is not surprising, since the 
l—— ', heaviest element in the (TMTSF) 2 X family is Se. Thus, unlike in the work of LMO, we consider here only the weak 

spin-orbit coupling limit corresponding to the equal spin triplet state (ESTP) proposed by one of us Ji7|-^9t. Here we 
, also take into consideration the effects of time reversal symmetry breaking into the structure of the order parameter, 

■ for finite magnetic fields. We would like to mention, in passing, that the ESTP state is also consistent with the 
OS \ absence of the Knight shift for fields parallel to the y (&') |l4| , predicting in fact that the low field susceptibility in the 

superconducting state at zero temperature Xs(T = 0) = Xn, the susceptibility at the normal state for any direction 
of the applied field 

We focus in this paper only on the angular dependence of the upper critical field, when the magnetic field is applied 
along the yz plane and show that deviations from perfect alignment along the y (&') axis cause a dramatic drop in 

■ the critical temperature for intermediate values of the magnetic field in the range of 5 to 30 Teslas. We then argue 
that this rapid suppression of the critical temperature is in qualitative agreement with the rapid drop of the putative 
upper critical in (TMTSF) 2 C10 4 @]. 

We model quasi-one-dimensional systems via the energy relation 



£ a ,<r(k) = e a (k) - afi B H, (1) 

with the a-branch dispersion 

e Q (k) = vp(ak x — kp) + t y cos(k y b) + t z cos(k z c), (2) 

corresponding to an orthorombic crystal with lattice constants a, b and c along the x, y and z axis respectively. 
In addition, since Ep = vpkp/2 t y t z , the Fermi surface of such systems is not simply connected, being open 
both in the xz plane and in the xy plane. Furthermore, the electronic motion can be classified as right-going (a = +) 
or left-going (a = — ). 

We begin our discussion by analysing the electronic motion for magnetic fields H = (0, H y ,H z ) in yz plane. We 
work in units where the Planck's constant and the speed of light are equal to one, i.e., h = c* = 1, and, throughout 
the paper, we use the gauge A = (H y z — H z y, 0, 0), where a and k x are still conserved quantities (good quantum 
numbers) while k y and k z are not. In this gauge, the semiclassical equations of motion become 

x — xo + avFt, (3) 
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y = Vo + a(t y b/u Cy ) cos(k y (Q)b + au Cy t), (4) 

z = zq — a(t z c/tu Cz ) cos(fc z (0)c — auj Cz t), (5) 

with u> Cy — vpGy and w Cz — vpG z , while G y — \e\H z b and G z — \e\H y c, where H z = H cos9 and H y = H sm9, where 
9 is the angle between H and z. The frequencies u Cy and lu Cz are independent and characterize the periodic motion 
along the y direction and z direction, respectively. 

This semiclassical analysis suggests a classification of three different regimes of confinement for the electronic 
motion: (a) a one-dimensional regime, where the semiclassical amplitudes of motion are confined along both the y 
and z direction, i.e., t z /u Cz < 1 or t y /u Cy -C 1 (provided that there is no magnetic breakdown, and that Ep ^> u> Cy 
and Ep 3> oj Cz ); (b) a two-dimensional regime, where the semiclassical amplitudes of motion are confined along the z 
direction but not along the y direction, i.e., t z /to Cz <C 1 or t y /cu Cy 3> 1; and (c) a three-dimensional regime, where the 
semiclassical amplitudes of motion are not confined in either y or z directions, i.e., t y /oj Cy 3> 1 or t z /oj Cz 3> 1. The 
confinement (localization) first occurs along the z direction when \t z c/uj Cz \ < \t y b/uj c | which implies | tan#| > t z /t y , 
i.e., when 8 > 6 C — arctan(t z /t y ). The localization first occurs along the y direction when | tan#| < t z /t y or 6 < 9 C . 
Hence in the limit where t z /t y <C 1 the localization happens first along the z direction except for a very narrow region 
of angles 8 < 9 C . For instance, taking t z /t y w 0.1, 9 C w 5.71° and the angular region where localization takes place 
first along the z direction is very wide (9 C < 9 < 90°), except for a few degrees from the z axis. 

Now we turn our attention to the quantum aspects of the problem. In the presence of the magnetic field H, the 
non-interacting Hamiltonian is 

H (k- eA) = e Q (k- eA) - a\i B H (6) 
in the gauge A = (H y z — H z y, 0, 0). The eigenfunctions of H (k — eA) are 

$ g „(r) = exp[ik x x]J N n ( — ) Jn z -u z ( — ) , (7) 
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where r = (x, y, z), y = n y b and z = n z c with associated quantum numbers qn = a, k x , N y , N z , a and eigenvalues 

e 9 n = £a,<T(k p ) + aN y uj Cy + aN z u Cz . (8) 

The function J p (u) is the Bessel function of integer order p and argument u, while now 

e a<(T {kp) = v F (ak x — fcj?) — afj, B H (9) 

is a ID dispersion. Notice that the eigenspectrum in Eq. (||) involves many magnetic subbands labeled by the 
quantum numbers N y and N z and that the eigenvalue e qn is invariant under the quantum number transformation 
(a, k x , N y , N z , o) — * (—a,—k x ,~Ny,—N z ,a), for the same spin state a. In addition, these magnetic subbands are 
spin-split into spin-up and spin-down bands. Thus, Cooper pairs can be easily formed in a ESTP pairing state, 
involving electrons with quantum numbers (a, k x , N y , N z , a) and (—a, — k x , —N y , —N z , a), provided that the pairing 
interaction conserves spin (which seems to be the case for the (TMTSF^X family, except for very small spin-orbit 
and dipolar couplings) p2fl . Thus, the order parameter vector 

A(r) = [A +1 (r),0,A_ 1 (r)] (10) 

has components in the m s = +1 (fx =TT) arL d m s = —1 (/i =[[) channels, only. 
In the case of ESTP, the gap equation can be written as 



A„(r) = A M A nInA^ + S n v n z ) (11) 



where r = (x, y, z), Sjv„iV e — (0, N y b, N z c). The matrix operator 

A %l*r. = E N y N x {x)FM Nm (q x - K NyNz ) (12) 
with Kn v n z = 2N y Gy + 2N Z G Z involves two terms, the prefactor 
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and the operator function 



E NyNz (x) = ex P [-^NyN^)] 



with Qmw — (Miy + M 2y )G y + {M\ z + M2 Z )G Z , and q x = —id/dx. The weight function 

W^{M iv ) = r% viN ^M 2u )T% v ^(M lv ), 
contains the single-particle renormalization factor 

T% v ,n,(M 2v ) = P% y (M 2y )P^(M 2z ) 

with coefficient 

P« (M 2v ) = J M2v (^) J M2 „-^ (— 

\ Ur„ / V UJ r „ 



(13) 
(14) 

(15) 
(16) 
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Notice that the second term in the right hand side of Eq. jl5|), i.e., N (Mi„) can be obtained from Eq. ( [l6| ) via 

the the following substitution a — > a and M 2v — > Mi,/. In addition, the operator Fn v nAQx) contains the Cooper 
singularity contribution 
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with Affi being the spin-dependent density of states, and Cl = (2cj ^7 /kT). 

The physical interpretation of Eq. (|ll]) is as follows. The order parameter A fJi (x, y, z) at a given position in space 
couples with the order parameter at A M (:c, y + N y b, z + N z c) at a different position, thus indicating that Josephson 
coupling between chains oriented along the x direction occurs. The strength of this coupling is controlled by the 
matrix operator A£ N , which is dependent both on the magnitude of the magnetic field and its direction, through 
G y and G z . This result is extremely interesting because it shows that the quantization effects of the magnetic field 
can make quasi-one-dimensional superconductors evolve from a strongly coupled Ginzburg-Landau local regime (low 
magnetic fields) into a variable range Josephson coupled non-local regime (high magnetic fields). 

We will now focus our attention to three limits, where the critical temperature as a function of the applied magnetic 
(or conversely the the upper critical field as a function of temperature) can be calculated analytically. We will assume 
that the interaction A M = A and the density of states Af^ — Af are independent of the spin channel P3l , which implies 
that A^ = A in the analysis that follows. 

One- dimensional regime: As seen in the semiclassical analysis of Eqs. (||), (0) and (||), there is double confinement 
of the orbital motion when t y /u> c <C 1 and t z /uj c „ <C 1. Using the fact that q y and q z are good quantum numbers 
one can write 



A(x, y, z) = exp[i(q y y + q z z)]u(x), 
which substituted in Eq. ( |Tl| ) transforms it into the generalized (two-term) Hill equation 

d 2 
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x dx 2 



Aid - ^2 B y cos0„(x) 



u(x) = 0, 



(18) 



(19) 



to leading non-trivial order in q x . Here 4>v{ x ) = {<lvSv — 2G v x), v = y,z, s y — b and s z — c. The coefficient of the first 
term is ftid = (1 — ty/u^ — t 2 /uo 2 J. This coefficient renormalizes the characteristic length scale along the x direction 
l x = \/Cvf/T, which contains the constant C — 7£(3)/167r 2 . The second coefficient is 



flu la (ft) +J2 
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while B v — (t v juj c J) In \ r yuj Cv /2wT\ are the last coefficients corresponding to the amplitude of co-sinusoidal spatial 
oscillations. Physically, Eq. (|l9|) corresponds to a limit of weakly coupled Josephson chains, with magnetic field 
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dependent Josephson couplings B v . The highest critical temperature in high magnetic fields, to lowest order in 
t„ I , occurs when q y = q z 



0. Thus, T c is determined by the condition A\ d = 0, leading to 

2 
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with prefactor T Cld = (2u>dPf/ir) exp(— 1/XAf) p4j , since both the amplitude B v and the period ir/Gy of cos(2G v x) 
are very small. When t y /uj c <C 1 and t z /u> Cli <C 1 there is a magnetic field induced double quantum confinement 
(localization) effect along both the y and z direction, for almost all 9. However, localization along the y direction is 
harder to achieve in the Bechgaard salts (TMTSF^X (X = CIO4, PFg, ...), since t y is typically of the order of a 100K, 
thus requiring magnetic fields the order of 100T to take place. 

Two-dimensional regime: When the value of magnetic field is lowered and the angle 9 is changed to satisfy conditions 
t y /u! Cy ^> 1 and t z /uj Cz <C 1, we reach a two dimensional regime where quantum confinement occurs only along the 
z direction. In this two dimensional regime the effects of the magnetic field along the y direction can be treated 
semiclassically. Using Eq. (f^) the resulting eigenvalue equation is 



dx 2 



+ L y (x) + A 2d + B z cos <f) z (x) 



u(x) = 0, 



(21) 



where the coefficient of the first term is /?2d = (1 — ^z/^c )• The second term is L y (x) = (£ y (f> y (x)/b) 2 , with coefficient 
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Ctyb/Ty/2, while the third term is 

A 2 d = (3 2d HQ) 

The critical temperature resulting from Eq. ( pl| ) is 
1.1 
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where 
90.1 
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. T C =T, 

\i shown in Fig. (Q). 
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plot oo : T c as a function of magnetic field for various angles 9 close to 
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(22) 
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FIG. 1. This figure shows that the suppression of the critical temperature T c for applied fields H along the yz plane in the 
two-dimensional regime. The direction 6 = 90.0° corresponds to perfect alignment along the z axis. The parameters used were 
T C2d = 1.5K, t y — 100K, t z = 10K, with lattice constants characteristic of Bechgaard salts. 



Notice in Fig. (|l|), the rapid reduction in T c for very small deviations from 9 = 90.0° (H || y). A similar drop in 
T c has been seen experimentally in the putative critical temperature of (TMTSF) 2 C104 at high magnetic fields for 
small angular deviations from the y axis (b' axis) Jig ], 

Three-dimensional regime: The three dimensional regime is characterized by the absence of quantum confinement, 
and occurs only at low magnetic fields where t y /u> Cy ^> 1 and t z /uj Cz ^> 1. In this case the gap equation reduces to 
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u{x) = 



(23) 



where £ x and L y (x) were already defined, and L z (x) = (£ z (j) z (x)/c) 2 , with characteristic length scale along the z 
direction t z = \fCt z cjT\[2. The coefficient A^d = [in (Q) — l/AA/"] . This low magnetic field regime corresponds to 
the anisotropic Ginzburg-Landau limit, and thus the critical temperature is 



T c = T c (0) 



CV2 



(24) 



The angular dependence of the upper critical field in the xy plane (instead of yz plane discussed here) was studied 
by Huang and Maki in the Ginzburg-Landau regime. 

In summary, we have assumed an ESTP state as a plausible candidate for triplet superconductivity in quasi-one- 
dimensional systems [[17-19 , and we have presented analytical results for the angular dependence of the upper critical 
field (critical temperature) of these superconductors when the magnetic field is applied in the yz plane. We have 
discussed three general regimes: (a) a one-dimensional regime at very high magnetic fields where there is quantum 
confinement both along the y and z directions (t y /uj Cy -C 1 and t z /u) Cz <C 1); (b) a two-dimensional regime where 
there is quantum confinement only along the z direction (t y /tu Cy 3> 1 and t z /oj Cz < 1); and (c) a three-dimensional 
regime at low magnetic fields (Ginzburg-Landau limit, where there is no quantum confinement (t y /u) c S> 1 and 
t z /uj Cz ^S> 1). The one-dimensional limit is reached for quite high magnetic fields H > 100T, which are not available 
today from continuous field sources. However, the two-dimensional regime, which is in the range of 5 to 30T, is 
attainable with continuous field sources. This two-dimensional regime is very interesting from the experimental point 
of view, because there is a very rapid suppression of the predicted reentrant superconducting phase (for H || y) [|l]~[|, 
when the applied magnetic field is just fractions of a degree away from perfect alignment with the y (&') axis. Our 
analytical calculations seem to be in qualitative agreement with the possible experimental observation of such an 
effect in (TMTSF)2C104 ]l5|] . We also have shown that the three dimensional regime just corresponds to the usual 
anisotropic Ginzburg-Landau limit. We would like to thank the Georgia Institute of Technology, NSF (Grant No. 
DMR-9803111) and NATO (Grant No. CRG-972261) for financial support. 
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